Viscoelasticity and metastability limit in supercooled liquids 
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A supercooled liquid is said to have a kinetic spinodal if a temperature T sp exists below which 
the liquid relaxation time exceeds the crystal nucleation time. We revisit classical nucleation theory 
taking into account the viscoelastic response of the liquid to the formation of crystal nuclei and 
find that the kinetic spinodal is strongly influenced by elastic effects. We introduce a dimensionless 
parameter A, which is essentially the ratio between the infinite frequency shear modulus and the 
enthalpy of fusion of the crystal. In systems where A is larger than a critical value A c the metastability 
limit is totally suppressed, independently of the surface tension. On the other hand, if A < A c a 
kinetic spinodal is present and the time needed to experimentally observe it scales as exp[w/(A c — A) 2 ], 
where ui is roughly the ratio between surface tension and enthalpy of fusion. 
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When a liquid is cooled below its freezing point with- 
out forming a crystal, it enters a metastable equilibrium 
phase known as supercooled 0. A supercooled liquid 
is squeezed in an uncomfortable time region: if we are 
too fast in measuring its properties, the system cannot 
thermalize and an off-equilibrium glass is formed; on the 
other hand, if we are too slow, the system has the time 
to nucleate the solid, and we obtain an off-equilibrium 
polycrystaL and, eventually, a thermodynamically stable 
crystal 0, Q • What is the maximum degree of supercool- 
ing a metastable equilibrium liquid can reach ? 

The tricky point about this question is that it mixes 
aspects of the experimental protocol, with intrinsic prop- 
erties of the system. If we stick to cooling the system 
linearly in time, there is a minimum cooling velocity be- 
low which the system is bound to crystallize 3]. This 
velocity is inversely proportional to the minimum nucle- 
ation time as a function of temperature U we cannot 
cool slower than this minimum cooling rate, otherwise 
crystallization occurs. On the other hand, as we cool, 
the relaxation time increases steeply, and therefore the 
system necessarily leaves the supercooled phase, and be- 
comes a glass, at the temperature where the relaxation 
time becomes too large for this minimum cooling rate. 

To penetrate deeper in the supercooled region, one can 
use an ad hoc nonlinear cooling protocol: cool fast close 
to the temperature where crystallization is a concern Q , 
and where relaxation time is still small; and slow down 
at lower temperatures, to cope with the increasing relax- 
ation time, once the nucleation time starts raising again. 
Therefore, we may think that the unique limitation to the 
extent of supercooling is given by our capability of cool- 
ing slow and fast enough a sample, and that in principle 
there is no bound to supercooling a metastable equilib- 
rium liquid. 

In fact, our experimental capability is not the only li- 
mitation to supercooling a system. If at a certain tem- 
perature the relaxation time of the liquid tr exceeds the 
nucleation time of the crystal tn, no equilibrium mea- 



surements can be performed on the liquid sample and the 
supercooled phase does not exist anymore. Such a tem- 
perature is called kinetic spinodal T sp , and it marks the 
metastability limit of the supercooled phase 0, 0] ■ Such 
a metastability limit does not depend on the cooling pro- 
tocol, and is rather an intrinsic property of the sample. 
Below T sp the only equilibrium phase (either stable or 
metastable) is the crystal, whereas glassy and polycrys- 
talline off-equilibrium configurations can be obtained if 
we cool fast or slow enough, respectively. The aim of 
this work is to answer the following question: what is the 
main mechanism determining whether or not a metasta- 
bility limit is present in a given liquid ? 

Our answer to the question above is that this mecha- 
nism is viscoelasticity. The central idea is that on time 
scales shorter than the structural relaxation time, deeply 
supercooled liquids exhibit a solid-like response to strains 
In particular, the inclusion of a crystalline nucleus in 
the liquid matrix produces a strain 0,0, and thus a long- 
range elastic response 0, EJ EH E3 • Therefore, the ther- 
modynamic drive for crystal nucleation gets depressed by 
an elastic contribution, which depends on the ratio be- 
tween the relaxation time tr and the time scale tn over 
which nucleation occurs. In this way a self-consistent 
equation involving nucleation and relaxation times is ob- 
tained. 

We will show that liquids can be classified according to 
the magnitude of their elastic response relative to their 
bare thermodynamic drive to crystallization (that is the 
liquid-crystal Gibbs free energy density difference SG). 
We find that elastic effects can be large enough to sup- 
press the metastability limit, opening the possibility to 
supercool a liquid down to the point where its entropy 
reaches the entropy of the solid, the so-called Kauzmann 
temperature 7k [2|. Moreover, even if a metastability 
limit is present, the viscoelastic response determines how 
large is the relaxation time at the kinetic spinodal, and 
it thus establishes whether the metastability limit is ac- 
cessible within the maximum experimental time, or is 



2 



hidden by the off-equilibrium glassy phase. 

According to classical nucleation theory (CNT) the 
Gibbs free energy barrier to nucleation is given by 



is due to the formation of a crystal nucleus, and this 
happens on a time scale of the order of the nucleation 
time. Therefore we must set t — tn, and write, 
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where d is the dimension, a is the surface tension, 5G is 
the (positive) Gibbs free energy difference per unit vol- 
ume between the two phases, and a is a numerical con- 
stant (for a spherical nucleus in d = 3, a = 167r/3). 

When the nucleus of the stable phase forms in an elas- 
tic background, there is an extra energetic price that has 
to be paid Because of the long-range nature 

of the elastic strains, the elastic price is proportional to 
the nucleus volume, and thus it corrects the bare ther- 
modynamic drive 5G, 
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In a perfect isotropic solid, assuming that the elastic pa- 
rameters are the same in the two phases, the elastic en- 
ergy density is given by £ c i as tic = £o G m where G^ is 
the (constant) elastic shear modulus, and where [H|. 
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is a dimensionless constant: K is the bulk modulus and 
Sv/v is the stress-free volume misfit between the two 
phases. In a viscoelastic liquid we have a stress-relaxation 
function, or time-dependent shear modulus, G(t) [lflj . In 
this case, Schmelzer and coworkers have shown that the 
time-dependent elastic contribution can be written as, 
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where the infinite frequency shear modulus is, _G C 
G(t = 0), and the dimensionless function / is 
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f(t) = - I dt' G(t')/G c 
1 Jo 



(5) 



For a Maxwell liquid [3 G(t) = G^e^/^ and thus, 
f(t/T R ) = (1 - e- t/TR )T R /t. We shall not assume the 
Maxwell form, but we shall still make the hypothesis 
that the stress-relaxation function G depends on the ratio 
V T R, where tr is the liquid's structural relaxation time. 
In this case also / = /(*/tr). In log scale the function 
/ typically has a step-like shape, going from 1 = /(0) 
to = /(oo) with a sharp drop at /(l). The physical 
meaning is clear: for times much larger than the relax- 
ation time the liquid is able to relax the stress, and there 
is no elastic contribution, while for times much smaller 
than tr it responds as a solid, with finite shear modulus. 

What is the time t we have to plug into Elastic (^/ t r) ? 
The viscoelastic contribution to the nucleation barrier 



tn = tn ex P 



k B T [SG(T) - Elastic (t N /t R )] 
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From this formula we clearly see that the main effect 
of the elastic contribution is to increase the nucleation 
time. This effect becomes dramatic in solid-solid phase 
transitions, where nucleation may be totally suppressed, 
even with a nonzero thermodynamic drive SG, whenever 
SG(T) < £ clastic = eoGoo HSGSGII. In liquids, how- 
ever, the stress is always relaxed for long enough times, 
and eq.© always admits a solution tn(T), as long as 
tr < oo. 

In order to solve eq. ©, we have to discuss the In- 
dependence of the various quantities. For the relaxation 
time we shall assume a Vogel-Fulcher-Tamman (VFT) 
form, which is known to describe very well a wide variety 
of supercooled liquids on a large range of times Q , 
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This form of the relaxation time implies a singularity at 
7k, which is normally identified with the Adam-Gibbs 
thermodynamic transition and it is approximately 
equal to the temperature where the liquid-crystal entropy 
difference SS vanishes. The Gibbs free energy difference 
is often fitted to a linear form, 8G = (1 — T/T m )6h/v, 
where T m is the melting temperature, Sh is the molar en- 
thalpy of fusion, and v is the molar volume of the crystal 
[3 Ha. However, recalling that SG = -d{5S)/dT, con- 
sistency with |(7J requires that SG must reach its maxi- 
mum at Tk. Therefore, we shall write, 



SG(T) = g(T/T m ) (1 - T k /T m ) 8h/v 



(8) 



where the dimensionless function g satisfies the relations: 
g(l) = and g(Tk/T m ) = 1. The surface tension a will 
be assumed to be a constant, which is a reasonable hy- 
pothesis within CNT [T3|. 

The prefactor of nucleation 7$ in © is significantly 
T-dependent. According to CNT, is proportional 
to the product of an elementary thermal time scale, 
ftpiank/^B? 1 , and an Arrhenius term due to the barrier 
B a particle has to overcome to cross the nucleus in- 
terface, exp( J B/fc B T) [13 . It is customary to approxi- 
mate such term by 1/D (D is the self-diffusion coeffi- 
cient), and by using the Stokes-Einstein (SE) relation, to 



write r, 



1/D 



tr 



However, as the system be- 



comes viscous, the SE relation breaks down |18t HH Eof. 
and typically 1/D <C tr. Indeed, in the extreme case 
of solids, even though the shear viscosity (and thus tr) 
is infinite, D remains nonzero. Thus, the prefactor t^ 
remains much smaller than the relaxation time at low 
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temperatures, and it cannot by itself grant the fact that 
tn > m- 

A second delicate point about t n is its dependence on 
the sample's volume. The quantity studied by CNT is 
not directly the nucleation time, but the nucleation rate 
J, that is the number of nuclei per unit time, per unit 
volume. J is a constant for large enough samples and 
therefore the nucleation time of a given sample of volume 
V scales as 1/V 3]. Taking very small samples is indeed 
a well known trick to increase nucleation time. However, 
there is a limit to this procedure, given by the volume vq 
below which the 1/V scaling breaks down: below vq the 
nucleation rate is not anymore constant, either because 
surface effects become dominant over bulk properties, or 
because the sample size is smaller than the critical nu- 
cleus j2l| • The theoretical metastability limit of a liquid 
can be defined as the point where the largest possible 
nucleation time is surpassed by the relaxation time, and 
therefore we set tn = 1/(vqJ). Anyhow, we recall that 
the results are rather insensitive to the choice of the ref- 
erence volume, since this enters only in the prefactor r N , 
and not in the most relevant exponential factor. 

We measure all times in equation (j^J in units of r N , 
and define x = log(TN/T N ) and y = log(rpt/r^)). In this 
way equation JBJ can be rewritten as, 



T'x = lu [g{T')-Xf{e x /ey) 



where T' = T/T m is the reduced temperature, and, 
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Two more dimensionless parameters are present in the 
theory: they are the rescaled Kauzmann temperature 
T{ = T k /T m and the rescaled VFT barrier A' = A/T m . 
They both enter in the relaxation time, and thus in 
y = y(T'/Tl, A'). If we set A = 0, T k = and 
g = (1 — T/T m ), we recover the classic expression for 
the nucleation time, where exp(w) gives the scale of the 
minimum nucleation time (see, for example, H3)- O n 
the other hand, for A ^ we have elastic corrections to 
the nucleation time. 

Equation can be easily studied numerically once the 
forms of g and / are specified. However, the key physical 
features can be worked out in full generality. First, we 
note that for T — > T m , g — > whereas y remains finite, 
and thus the solution x of © diverges. As expected, 
tn — > oo at the melting point. In the opposite limit 
T — > Tk a graphical study of equation © shows that for 
A < 1 the nucleation time remains finite, while y — > oo, 
and thus a metastability limit T sp > T k exists. On the 
other hand, for A > 1, tn — > oo for T — > T k : in this 
case both relaxation and nucleation times diverge at T k , 
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FIG. 1: Nucleation and relaxation vs reduced temperature 
T = T/Tm, for different values of the parameter A. Full 
line, nucleation: x(T') = Ic^tn/tn)- Broken line, relaxation: 
y(T) — log(TR./Tf})). LQ stands for equilibrium liquid, GL for 
off-equilibrium glass, and CR for off-equilibrium polycrystal. 
Left panel: for A = 2, larger than A c = 1.58, there is no kinetic 
spinodal. Right panel: for A = 0.3, smaller than A c , a kinetic 
spinodal exists at T sp . Inset: for A = 0.03 the difference 
between glass and polycrystal is particularly ill-defined. The 
other parameters are: uj = 4, A/T m = 4, Tk/T m = 0.2. For 
f(t/Tn) we have taken a Maxwell form, while g(T') is a cubic. 



but they still may cross at a kinetic spinodal. To check 
this last possibility we set x = y in equation (|5J), and see 
whether there is a solution T sp > T k of this equation, 



T Q 



A' 



V - Tl 



wbcn-A/a)] 1 - 
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In this formula the factor SI = log(rp/r^) is not relevant, 
since it is large in modulus only at low temperatures, 
where however the factor (T" — T k ) _1 dominates. We 
recall that, by construction, /(l) is a number close to 
1/2 (in Maxwell's case /(l) = 1 - 1/e = 0.63). The 
l.h.s. of equation (|12|) diverges at T£, whereas the r.h.s. 
is finite at and diverges at a temperature larger than 
T{, defined by g(T') - A/(l) = 0. Given that g(TQ = 1, 
we conclude that for A < 1//(1) there is a solution T s ' p , 
and thus a metastability limit, while for A > 1//(1) there 
is no solution in the physical interval g(T') — A/(l) > 0. 

We recall that A is, up to factors of order one, the ra- 
tio between the characteristic scales, Goo and Sh, of the 
two competing mechanisms at work: the elastic forces, 
depressing nucleation, and the thermodynamic drive to 
crystallization, enhancing nucleation. What we have 
found is a critical value, 



A c = 1//(1) 



(13) 



separating systems with a kinetic spinodal, from systems 
without a kinetic spinodal. In those liquids where A < A c 
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the metastability limit may be shifted at lower T by elas- 
ticity, but it still exists, while for A > A c elastic effects 
are strong, and the nucleation time is amplified enough 
to completely suppress the kinetic spinodal. This dis- 
tinction is, of course, only valid to logarithmic accuracy: 
whenever A ~ A c , the nucleation time may be only slightly 
larger than the relaxation time, and it is wise to assume 
that the liquid is not really stable against crystallization. 

In Fig.l we show solutions of cq.® for A > A c (left 
panel), and A < A c (right panel). In the first case, there 
is no spinodal, and thus the three phases, glass, super- 
cooled liquid, and polycrystal, are well defined. In partic- 
ular, the off-equilibrium glassy and polycrystalline states 
are well separated from each other by the 'buffer' liquid 
phase. When A < A c , however, the situation is rather 
different. The equilibrium liquid phase is still well de- 
fined as long as it exists, that is for T > T sp ; however, 
it is hard to say what is the physical difference between 
the two off-equilibrium phases, i.e. glass and polycrystal. 
Although it may be operationally meaningful to distin- 
guish glass from polycrystal when T 3> T sp , there seems 
to be no obvious qualitative difference when T ~ T sp , or, 
even worse, when T < T sp . Close to the spinodal point, 
whether the system falls out of equilibrium by crossing 
the broken or the full line, is irrelevant. Glass and poly- 
crystal blend into a single off-equilibrium phase |22j. 

In systems with A > A c the supercooled phase is well 
defined down to Xk, where the relaxation time diverges. 
In this case, a solution of Kauzmann paradox (for exam- 
ple, in terms of a thermodynamic transition at Tk) must 
be found. On the other hand, when a spinodal is present 
the supercooled liquid ceases to exist at T sp , which is 
Kauzmann's resolution of the paradox 0|. 

When a kinetic spinodal exists, is it within the bound- 
aries of experimental observation ? To answer this ques- 
tion we have to define a spinodal time, r sp = r K (T sp ), 
and compare it to the maximum time of the experiment, 
r cxp . From eq. H12|) we obtain (for d = 3), 



' sp 



exp 



(1 - A/A c )2 



(14) 



Therefore, as expected, the position of the kinetic spin- 
odal depends also on the parameter u>; a large value of 
ui implies a large value of the surface tension a com- 
pared to the drive to crystallization 5h, and thus a larger 
nucleation time Therefore, when uj is significantly 
larger than (1 — A/A c ) 2 , the spinodal time may be too 
large compared to the experimental time, and thus the 
metastability limit be experimentally unaccessible |2.'"il| . 

Conventionally, a good glass-former is a system which 
does not crystallize easily under a linear cooling, i.e. a 
system whose minimum cooling rate r m i n is not too large. 
Given that r m j n ~ l/r m i n , where Tmin is the minimum nu- 
cleation time as a function of T j4j, the larger r m i n , the 
better the glass-former. To what extent T ln ; n depends on 
A ? When A <C A c the kinetic spinodal typically takes 



place at a point where the nucleation time is still a de- 
creasing function of T (Fig. 1, inset of right panel). In 
such case T min is effectively given by the spinodal point. 
From lfT4*jl we therefore have r m i n ~ exp[— w/(l — A/A c ) 2 ], 
and the glass-forming capability of the system is strongly 
dependent on A. On the other hand, when A ~ A c 
or A > A c , rmin is given as usual by the minimum of 
the nucleation curve. In this case T mm ~ exp(w), the 
glass- forming properties depend almost exclusively on lo, 
whereas their dependence on A is weak. 

In this work we studied the effects of viscoelastic re- 
sponse on the metastability limit of supercooled liquids. 
We defined a parameter A whose value rules whether a 
kinetic spinodal exists or not. This parameter encodes 
the competition between elastic response and thermody- 
namic drive to crystallization: for values of A below A c 
elasticity is relatively weak, the drive to crystallization 
always wins, and a metastability limit is present. For 
A > A c the elastic response is so large that nucleation 
in inhibited compared to relaxation, and a metastability 
limit is not present. 

We thank A. Angell, G. Biroli, C. Di Castro, I. Giar- 
dina, T.S. Grigera, and P. Verrocchio for many useful 
comments and clarifying discussions. 
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